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摘要 : 利用 载波 相位 偏差 产品 对 卫星 端 硬件 偏差 进行 改正 ， 采 用 星 间 单 差 消 除 接收 机 端 偏 差 ; 
结合 对 应 的 GPS 精密 轨道 和 钟 差 产 品 ， 依 次 固定 宽 埠 和 窄 埠 模糊 度 ， 然 后 得 到 基于 整数 解 的 
无 电离 层 组 合 模糊 度 ， 再 利用 虚拟 观测 值 对 非 差 无 电离 层 模糊 度 进 行 约束 ， 最 终 得 到 固定 解 轨 
道 。 同 时 ， 利 用 载波 相位 残 差 值 开展 星 载 天 线 相 位 中 心 变化 (PCV) 在 轨 估 计 ， 以 进一步 提升 定 
轨 精 度 。 分 别 利用 欧洲 定 轨 中 心 (CODE)、 法 国 空间 研究 中 心 (CNES) 和 武汉 大 学 (WHU) 发 
布 的 GPS 精密 轨道 、 钟 差 和 相位 偏差 产品 ， 对 GRACE-FO 卫星 进行 模糊 度 固 定 和 精密 定 轨 ， 
评估 不 同 机 构 产 品 的 性 能 。 结 果 表 明 : 天 线 PCV 改正 能 够 有 效 减 小 载波 相位 残 差 ， 提 升 定 轨 
精度 ; 采用 不 同 机 构 产 品 的 GRACE-FO 卫星 宽 埠 模糊 度 固定 率 均 优 于 99%， 窜 奉 模 糊 度 固定 
率 优 于 9596; 对 基于 CODE, CNES fll WHU 产品 的 GRACE-FO 卫星 精密 定 轨 结果 进行 评估 ， 
与 JPL 精 密 轨 道 相 比 ， 固 定 解 轨道 三 轴 平 均 精 度 优 于 7.0 mm; SLR 测 距 残 差 RMS 值 分 别 优 于 
9.6, 10.7 和 9.1 mm; 与 KBR 星 间 测 距 值 相 比 ， 利 用 固定 解 轨 道 计算 的 星 间 距离 精度 分 别 优 于 
1.8, 2.3 和 2.1 mm。 利 用 不 同 机 构 产 品 解 算 的 GRACE-FO 卫星 轨道 整体 差异 小 于 2 mm. 
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1 引 言 


低 轨 卫星 及 其 编队 在 遥感 、 海 洋 环境 监测 、 时 变 重力 场 估计 、 地 球 磁场 研究 、 掩 星 
气 观测 、 定 位 导航 授时 等 一 系列 科学 研究 和 工程 应 用 中 发 挥 着 不 可 或 缺 的 作用 。 低 轨 卫 
星 绝 对 精密 轨道 和 相对 位 置信 息 是 完成 上 述 任 务 的 前 提 和 基础 。 自 全 球 定位 系统 (global 


positioning system, GPS) 成 功 应 用 于 TOPEX/POSEIDON 卫星 精密 定 轨 后 "”， 基 于 星 载 
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GPS 观测 数据 的 简化 动力 学 定 轨 技术 "已 广泛 应 用 于 重力 场 恢 复 和 气候 试验 卫星 "(Gravity 
Recovery and Climate Experiment, GRACE) 、 挑 战 性 微型 卫星 载荷 ””(Challenging Mini- 
satellite Payload, CHAMP)、 海 洋 ”、 风 云 ”等 卫星 的 轨道 确定 。 利 用 GPS 精密 轨道 和 钟 
差 产 品 ， 低 轨 卫 星 模糊 度 浮 点 解 定 轨 精 度 可 达 1 ~ 3 cm  ””， 固 定 载波 相位 模糊 度 、 开 
展 接 收 机 天 线 相 位 中 心 变化 (phase center variation, PCV) 在 轨 估 计 能 够 进一步 提升 定 轨 

Bertiger 等 人 "和 Tapley A” AHER GPS 观测 数据 进行 TOPEX/POSEIDON E 
星 定 轨 ， 径 向 精度 达到 4 cm。 此 后 ， 多 个 低 轨 卫星 上 搭载 GPS 接收 机 用 于 精密 轨道 确 
Æ. Van den IJssel 等 人 ”和 赵 齐 乐 等 人 利用 国际 全 球 导 航 卫 星系 统 服务 (International 
Global Navigation Satellite System Service, IGS) 发 布 的 GPS 精密 星 历 ， 采 用 无 电离 层 组 
合 观测 值 进行 CHAMP 卫星 精密 定 轨 ， 径 向 精度 优 于 4 cm; Kang 等 人 "利用 卫星 激光 测 
EE (satellite laser ranging, SLR) 数据 对 GRACE 定 轨 精度 进行 评估 ， 三 维 轨 道 精度 为 2 和 ~ 
3 cm; Luthcke 等 人 "综合 利用 星 载 多 普 勒 无 线 电 定 轨 定 位 系统 (Doppler orbitography by 
radiopositioning integrated on satellite, DORIS). SLR 和 海面 测 高 数据 对 JASON-1 卫星 
GPS 简化 动力 学 定 轨 精 度 进行 评估 ， 满 足 1 cm 的 径 向 精度 要 求 ， 秦 建 等 人 "利用 星 载 观 
测 数 据 获 取 了 厘米 级 的 JASON-2 卫星 轨道 ， 在 其 他 低 轨 卫星 任务 中 ， 如 地 球 重 力 场 和 海 
洋 环流 探测 卫星 和 (Gravity Field and Steady-state Ocean Circulation Explorer, GOCE), 
TerraSAR-X“ 、 海 洋 ”、 风 云 ” 等 ， 也 得 到 了 类 似 结果 。 上 述 定 轨 中 模糊 度 参数 为 浮 点 
解 ， 为 实现 模糊 度 固定 ，Kroes 等 人 " "通过 构建 相对 动力 学 模型 并 固定 双 差 模糊 度 进行 
GRACE 双星 相对 定 轨 ， 再 利用 K 波段 测 距 (K-band ranging, KBR) 测量 出 星 间 距离 ， 发 
现 其 精度 优 于 1 mm. Jággi 等 人 "研究 了 固定 星 间 、 星 地 双 差 模糊 度 对 GRACE 定 轨 精度 
的 影响 ， 相 对 于 浮 点 解 ， 固 定 模糊 度 可 以 将 双星 相对 定 轨 精度 由 厘米 级 提升 至 1 mm; Zhao 
和 Hu “通过 固定 GRACE 双 差 模糊 度 实现 了 2 mm 的 相对 定 轨 精 度 。 但 是 对 于 单 颗 低 轨 卫 
星 ， 无 法 组 成 双 差 载波 观测 值 ， 限 制 了 该 方法 的 使 用 。 

为 实现 单 接收 机 模糊 度 固定 ， 国 内 外 学 者 先后 提出 了 小 数 周 偏差 法 "、 整 数 相位 钟 
法 ”和 人 解 厢 钟 差 法 ”。Bertiger 等 人 ”利用 美国 喷气 推进 实验 室 (Jet Propulsion Laborato- 
ry, JPL) 提供 的 载波 相位 偏差 产品 实现 了 单 星 模糊 度 固定 ，GRACE 相对 定 轨 精度 由 浮 点 
ES 解 的 6 mm 提升 至 2 mm, JASON-2 卫星 径 向 定 轨 精 度 优 于 1 cm. Montenbruck 等 人 ”在 
本 Sentinel-3A 卫星 定 轨 时 利用 法 国 国 家 空间 研究 中 心 (Centre National d'Etudes Spatiales, 

CNES) 发 布 的 宽 埠 改正 数 和 钟 差 产品 “对 模糊 度 进行 固定 ， 国 定 解 轨道 重 倒 弧 段 精度 和 
SLR 测 距 残 差 精 度 均 有 明显 提升 ， 张 小 红 等 人 ”基于 运动 学 方法 实现 了 GRACE 卫星 精密 
单 点 定位 (precise point positioning, PPP) 模糊 度 固 定 ，Arnold 等 人 "利用 欧洲 定 轨 中 心 
(Center for Orbit Determination in Europe, CODE) 发 布 的 相位 偏差 产品 实现 GRACE 和 
Sentinel-3 卫星 精密 定 轨 中 模糊 度 参 数 的 固定 ，Guo AP 对 比 了 双 差 模糊 度 固 定 、 单 星 模 
糊 度 固定 和 单 星 模糊 度 固定 基础 上 的 双 差 模糊 度 固 定 对 GRACE 卫星 定 轨 精度 的 影响 。 
另 一 个 影响 低 轨 卫星 定 轨 精 度 的 因素 是 星 载 接收 机 天 线 相位 中 心 的 位 置 。 低 轨 卫 星 发 
射 前 会 在 地 面 对 其 天 线 相 位 中 心 进行 测量 ， 但 是 由 于 在 轨 环 境 与 地 面 环 境 的 不 同 (如 多 径 ) 
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会 导致 天 线 相位 中 心 实际 值 与 标定 值 不 一 致 。Haines A” 利用 JASON-1 星 载 观测 数据 
计算 天 线 PCV， 经 PCV 改正 后 无 电离 层 组 合 载波 相位 残 差 由 8 mm 降低 至 5 mm， 基 于 定 
轨 结 果 计算 的 SLR 测 距 残 差 精度 优 于 11.9 mm. Jággi 等 人 ”对 GRACE 卫星 天 线 相位 中 
心 进行 了 研究 ， 并 描述 了 两 种 天 线 PCV 的 估计 方法 : 直接 估计 法 和 残 差 法 。Montenbruck 
55 AP" X Sentinel-3A 卫星 的 天 线 PCV 进行 了 在 轨 估 计 ， 国 内 马 洋 等 人 "对 JASON-2 卫 
星 的 天 线 PCV 进行 了 估计 ， 改 正 后 载波 残 差 减 小 至 5.5 mm， 径 向 定 轨 精度 优 于 1 cm. 
GRACE 后 续 星 (GRACE follow on, GRACE-FO) 卫星 于 2018 年 5 月 成 功 发 射 。 与 
GRACE 卫星 类 似 ，GRACE-FO 卫星 编队 由 两 颗 卫 星 组 成 (GRACE-C 和 GRACE-D P 
星 )， 星 上 均 搭载 了 GPS 接收 机 、SLR 角 反 射 器 和 KBR 设备 ”。Kang 等 人 ”研究 了 仅 
使 用 GPS 数据 和 联合 使 用 GPS 和 KBR 数据 对 GRACE-FO 定 轨 精 度 的 影响 ， 邵 凯 等 
AS 开展 了 基于 模糊 度 固定 的 GRACE-FO 精密 定 轨 研究 。 本 文 利用 星 载 实测 数据 开展 
GRACE-FO RA PCV 在 轨 建 模 ， 同 时 分 别 使 用 CODE, CNES 和 武汉 大 学 (WHU) 的 GPS 
精密 轨道 、 钟 差 和 相位 偏差 产品 进行 GRACE-FO 单 星 模糊 度 固定 及 精密 定 轨 ， 研 究 模糊 度 
固定 和 天 线 PCV 在 轨 标 定 对 低 轨 卫星 绝对 、 相 对 定 轨 精 度 的 影响 。 
第 2 章 在 GPS 观测 方程 的 基础 上 ， 给 出 单 星 模糊 度 固 定 和 天 线 PCV 估计 的 数学 模型 。 
第 3 章 介 绍 了 GRACE-FO 卫星 定 轨 策略 。 第 4 章 对 结果 进行 分 析 ， 包 括 天 线 PCV 计算 结 
果 、 模 糊 度 固定 效果 、 单 星 绝对 定 轨 精度 以 及 双星 相对 定 轨 精 度 。 第 5 章 给 出 了 结论 。 


2 ”数学 模型 


2.1 GPS 观测 模型 
GPS 伪 距 和 载波 相位 观测 值 可 表示 为 : 


nj 


c(d — dt*) — IS, + A(N2; + Bu; — BS) + Ago? + Apoo + Apcv 


| Ps = pi + c(dt, — dt’) + I2; + bej — bS + Apco + Apcv 


AF, Ps, LI. 表示 伪 距 和 载波 相位 观测 值 ，r 和 s 分 别 表示 星 载 接收 机 和 卫星 ，; 表示 
信号 频率 ，ps 为 接收 机 和 卫星 间 的 几何 距离 ，c Ati, dt, 和 dts 为 接收 机 和 卫星 钟 误差 ， 
I, 为 电离 层 延迟 ， 广 和 bs 为 接收 机 和 卫星 端 伪 距 硬件 偏差 ，B,; 和 Bs 为 载波 相位 硬件 
偏差 ， 以 周 表 示 。 Apco 和 Apov 表示 星 载 天 线 相位 中 心 偏差 (phase center offset, PCO) 和 
PCV 对 观测 值 的 影响 ， 其 中 Apcv 随 高 度 角 和 方位 角 变 化 。 入 ; 为 频率 7 对 应 的 载波 波长 ， 
NS, 为 相位 模糊 度 ，ws 为 载波 相位 缠绕 。 相位 缠绕 可 利用 模型 进行 改正 ”， 后 续 推 导 中 将 
忽略 该 项 。 

为 消除 电离 层 延 迟 一 阶 项 ， 须 构建 无 电离 层 组 合 观测 值 。 同 时 在 IGS 等 数据 处 理 中 ， 
通常 以 伪 距 观测 值 为 基准 解 算 导 航 卫 星 钟 差 ， 此 时 接收 机 和 卫星 端 伪 距 硬件 偏差 会 被 吸收 
进 钟 差 项 ， 式 (1) 可 写 为 : 


A 
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Phy = pr c(dt, — dt*) + Apco + Apcv 
L i = pf + cldt, — d£) + ANSqp+ Arco + Apcv — (2) 
Nsp = = Nis + (Bye — brIF/A) — (Bip — big/ X) 
AF, Pi 和 Lig 为 无 电离 层 组 合 伪 距 和 载波 相位 观测 值 ，dt. 和 d 为 吸收 了 伪 距 硬件 
延迟 的 接收 机 和 卫星 钟 差 ，Ns1j 为 无 电离 层 组 合 模糊 度 ， 和 Dre 为 接收 机 和 卫星 端 无 
电离 层 组 合 伪 距 硬件 偏差 ， Bi 和 Bip 为 组 合 载 波 相 位 偏差 ， 入 为 L 载波 波长 ，Nsr 为 
利用 伪 距 和 载波 相位 观测 值 进行 定位 时 解 算 的 浮 点 解 模 籼 度 ， 其 包含 了 接收 机 、 卫 星 端 伪 距 
和 载波 硬件 偏差 。 
2.2 ” 单 星 模 糊 度 固定 
无 电离 层 组 合 模 糊 度 NS ip n] 4HRUI e RU e CE RET. 


Vs fife s fi s 
Nip = B- gw j 万 十 万 NP wt 
Newt = Nh 一 人 rs ? 3) 


Nur = NEn t+ dur — di, 
HP, fi 和 fo 为 对 应 载波 L Al Lo 的 频率 ，Nswr 和 Nsw 为 整数 宽 埠 和 窄 埠 模糊 度 ， 
NesNr 为 浮 点 窄 埠 模糊 度 ， N31 和 Nea 为 载 ; X Li 和 了 2 的 模糊 度 ， dr NL 和 dxr, 为 接收 机 和 
卫星 端 伪 距 和 载波 硬件 偏差 的 窗 巷 组 合 。 利 用 MW HB T 解 算 宽 埠 模糊 度 ， 可 得 : 
a fia folo 万 Psi 一 户 Ps 

RP, dowi 和 dy, 为 接收 机 和 卫星 端 宽 埠 硬件 偏差 ，Awr 为 宽 埠 波长 。 分 别 对 接收 机 和 
卫星 端 宽 巷 偏差 进行 改正 ， 可 实现 宽 奉 模糊 度 固 定 。 宽 奉 模 糊 度 固 定 后 ， 代 入 式 (3) 中 计算 

FARREZ, TA: 
Ne 一 


fo fit fe 
fi Ue fo fi 
同样 ， 对 接收 机 和 卫星 端 罕 埠 偏 差 进行 改正 可 以 实现 鹤 巷 模糊 度 的 固定 。 实 际 处 理 时 ， 卫 星 
端 宽 埠 和 罕 埠 模糊 度 相 位 偏差 可 采用 地 面 站 网 数据 进行 解 算 并 向 用 户 提供 ， 用 户 端 宽 埠 和 
罕 巷 硬件 偏差 则 无 法 直接 获取 。 文 中 采用 星 间 单 差 的 方式 消除 接收 机 端 偏差 ， 根 据 式 (4) 和 
式 (5) 可 得 到 : 


New. Noir = Noni + drNL — dNL - (5) 


| ANY, = ANe Wi, — Adit? (6) 
ANNE = ANENE ~ Ady? 

UB. ANS M ANNE 为 卫星 st 和 so A. ARF AAPOR AY E Ae, ANS? 和 
ANS 为 整数 寅 、 FERE, Adv? 和 Ads? 为 星 间 宽 、 窄 替 硬 件 偏差 的 差 值 ， 
可 利用 改正 数 进 行 改正 。 "T EERME. IRIEN (7) 恢复 无 电离 层 组 
合 模糊 度 ， 即 : 


fife ANS ,82 fi 


png co Oe Fe fe 


(ANENE — Adni’) > (7) 
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AN 到 ”为 星 间 无 电离 层 组 合 模糊 度 的 差 值 。 利 用 式 (7) 对 非 差 无 电离 层 组 合 模糊 度 进行 约 
束 ， 得 到 固定 解 轨 道 。 
2.3 “天线 相位 中 心 估计 

GPS 相位 观测 值 是 通过 测量 导航 卫星 信号 发 射 时 刻 发 射 天 线 相 对 于 接收 时 刻 接收 天 线 
各 自 的 瞬时 相位 中 心得 到 ， 但 是 天 线 瞬时 相位 中 心 随 高 度 角 、 方 位 角 和 信号 频率 的 变化 而 变 
化 。 瞬 时 天 线 相位 中 心 可 用 相位 中 心 的 平均 值 (PCO) 和 相对 于 该 平均 值 的 变化 (PCV) 进行 
描述 ， 如 式 (8) Pra”: 


Aó(a, z, fi) = Apco + Apcv(a, z, fi) ， (8) 


AF, a 为 导航 卫星 相对 于 接收 机 的 方位 角 ，z AREA, fi 为 信号 频率 ，i = 1,2 表示 信 
号 频率 数量 。 利 用 上 式 可 实现 瞬时 相位 中 心 到 天 线 参 考点 的 转换 ，Apco 和 Apcv 在 使 用 时 
低 轨 卫星 星 载 天 线 相位 中 心 可 采用 球 谐 函 数 或 分 段 线性 函数 进行 描述 ” ”。 本 文 利用 
> 分 段 线性 函数 对 低 轨 卫星 接收 天 线 PCV 进行 建 模 。 该 模型 假设 天 线 PCV 由 不 同 格 网 构成 ， 
对 天 线 方位 角 和 高 度 角 进行 等 分 得 到 格 网 ， 其 中 方位 角 范 围 为 0" ~ 360*， 高 度 角 范围 为 
0° ~ 90"。 当 某 时 刻 观 测 值 处 于 格 网 ABCD 内 P 点 时 ， 则 其 PCV 值 采用 线性 内 插 为 ”: 


Apcv,p =(1 — a)(1— B)Apcv,A + o(1 — B)Apcv,a + aBApcv,o + (1 — a)8Apcv,p = 
[(1— o)(1 — 8) a(1— 8) a8 (1— o)&][Aecv,A Apcvp Apcvc Apcvpl (9) 


其 中 a= (a -— a)/(a — ai), 8 = (z — z)/(z2 — 21) Apov.p NP AN PCV fü, o $8 
为 组 合 系数 ，a 和 zz 为 PP 点 的 方位 角 和 高 度 角 ，ai NA AD ABI. a AB ^C 
点 的 方位 角 ， 2 为 4 和 B AREA, 2 ACAD AREA, Apcv.a, Apcv,B, Apcv c， 
Apcv,p 分 别 为 A, B, C, D 点 的 PCV 值 ， 为 待 估 参 数 。 文 中 利用 载波 相位 残 差 解 算 各 格 网 
点 PCV 值 ， 定 轨 时 再 利用 模型 内 插 出 特定 方位 角 和 高 度 角 下 的 PCV 改正 数 。 


3 定 轨 策略 


采用 简化 动力 学 方法 进行 GRACE-FO 卫星 精密 定 轨 ， 估 计 参 数 包括 卫星 初始 位 置 和 速 
度 矢量 。 利 用 卫星 属性 信息 ” 进行 太阳 光 压 和 大 气 阻力 建 模 ， 采 用 阻力 温度 模型 94 (drag 
temperature model 94, DTM94°") 计算 大 气 密度 ， 定 轨 过 程 中 同时 估计 大 气 阻 力 和 太阳 光 
压 系数 。 为 了 补偿 模型 误差 ， 计 算 切 向 和 法 向 每 圈 1 次 (one cycle per resolution, CPR) 的 
周期 性 经 验 加 速度 ， 其 他 待 估 参 数 还 包括 接收 机 钟 差 和 模糊 度 。 星 载 GPS 天 线 PCO 根据 
公布 的 参数 进行 改正 ”，PCV 利用 载波 相位 残 差 进行 在 轨 估 计 ””。GPS 卫星 轨道 和 钟 差 
采用 事后 精密 产品 ， 同 时 利用 相位 偏差 产品 对 载波 相位 观测 值 进行 改正 。 定 轨 时 首先 基于 无 
电离 层 组 合 观测 值得 到 浮 点 解 轨道 ， 再 根据 第 2 章 介 绍 的 方法 ， 依 次 固定 宽 埠 和 窒 巷 模糊 
度 并 恢复 无 电离 层 组 合 模 糊 度 ， 最 终 得 到 固定 解 轨道 。 具 体 定 轨 策略 见 表 1。 
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#1 GRACE-FO 卫星 定 轨 和 模糊 度 固 定 策略 


模型 参数 说 明 
地 球 重力 场 EIGEN.06C(130X 130) 
固体 潮 和 极 潮 IERS 2010 
Tj id FES2004(30 X 30) 
m 海洋 极 潮 Desai (30X 30 
先 验 力 模 型 三 体 引力 JPL e bo 
相对 论 效 应 IERS 2010 
太阳 光 压 [Lg ZUM 
大 气 阻 力 卫星 宏 模 型 ，DTM94 
观测 数据 非 差 无 电离 层 组 合 
采样 间隔 10s 
截止 高 度 角 0? 
GPS 轨道 CODE, CNES, WHU 产品 
使 用 数据 和 产品 “GPS 钟 差 和 相位 偏差 CODE, CNES, WHU 产品 
GPS 天 线 igs14.atx 
接收 机 天 线 改正 
相位 缠绕 改正 
相对 论 改正 IERS 2010 
TERS 卫星 初始 位 置 和 速度 矢量 
大 气 阻 力 系数 每 90 min 计算 1 次 
太阳 光 压 系数 每 定 轨 弧 段 计算 1 个 
参数 估计 经 验 力 1-CPR 切 向 、 法 向 经 验 加 速度 ， 每 90 min 计算 1 组 
接收 机 钟 差 每 历 元 计算 
模糊 度 每 跟踪 弧 段 计算 1 个 
—À Se p xE BUB 0.25 Jj. 7265/8] BUR 0.15 JH], 
人 高 度 角 阔 值 为 3 


4 ” 定 轨 结果 分 析 


GRACE-FO 卫星 轨道 高 度 为 500 km， 两 颗 卫 星相 距 约 200 km. 本文 使 用 2019 FŒ 
载 GPS 观测 数据 进行 GRACE-FO 卫星 精密 定 轨 。 观 测 数据 可 从 德国 波茨坦 地 球 科学 研究 
中 心 (Geo Forschungs Zentrum, GFZ) 提供 的 Level 1B RL04 数据 包 中 获取 。 数 据 包 内 同时 
含有 卫星 四 元 素 姿态 信息 、 机 动 信息 、JPL 提供 的 精密 轨道 以 及 KBR 星 间 测 距 值 等 。 文 中 
将 利用 IPL 精密 轨道 和 KBR 数据 进行 定 轨 精 度 评 估 。 分 别 利 用 CODE ^, CNES” ”和 
WHU ”发 布 的 GPS 卫星 精密 轨道 、 钟 差 和 相位 偏差 产品 进行 GRACE-FO 单 星 模糊 度 
定 和 精密 定 轨 。 

4.1 ”天线 相位 中 心 校正 

根据 2.3 节 描 述 ， 利 用 无 电离 层 组 合 载波 相位 残 差 计算 GRACE-FO 星 载 GPS RK 
PCV 改正 数 。 为 减 小 计算 量 ， 方 位 角 和 高 度 角 分 辨 率 均 为 5?。 具 体 计算 时 ， 首 先 利用 单 
个 定 轨 弧 段 的 残 差 数据 估计 该 弧 段 对 应 的 PCV 改正 数 ， 之 后 对 2019 年 所 有 有 效 的 PCV 
改正 数 进行 平均 ， 从 而 得 到 整个 定 轨 时 段 的 PCV 估计 值 。 对 载波 观测 值 进 行 PCV 改正 并 
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再 次 进行 定 轨 解 算 ， 利 用 新 的 残 差 数据 估计 PCV 改进 值 ， 共 迭代 4 次 ， 得 到 最 终 的 天 线 
PCV 值 。 图 1 为 GRACE-FO 卫星 的 天 线 PCV 计算 结果 ，a) 为 GRACE-C 卫星 的 PCV， 
b) 为 GRACE-D 卫星 的 结果 。 对 于 GRACE-C 卫星 ， 方 位 角 0? 为 卫星 飞行 方向 ， 对 于 
GRACE-D 卫星 ， 方 位 角 180° 为 飞行 方向 。GRACE-C 卫星 在 方位 角 109. AE f 509 Hj, 
PCV 有 最 大 值 13.28 mm; 在 方位 角 185° AEH 25° IN, PCV 存在 最 小 值 -12.57 mm. 
GRACE-D 卫星 在 方位 角 190"、 高 度 角 50° 时 ， 存 在 最 大 值 11.62 mm; 在 方位 角 5°、 高 度 
fg 25° It, PCV 存在 最 小 值 —12.15 mm. 


30° 15. 
10 

> > 

0 B © 

(am n 

-5 X R 

K K 
10 


注 : a) GRACE-C BÉ; b) GRACE-D 卫星 。 
1 GRACE-FO DA 5?X5? 天 线 PCV 改正 


利用 上 述 天 线 PCV 计算 结果 对 载波 相位 观测 值 进行 改正 ， 再 次 进行 定 轨 处 理 ， 对 载波 
相位 残 差 值 进行 统计 ， 得 到 残 差 均值 如 图 2 所 示 ， 此 时 绘图 阀 值 为 5 mm ARR PCV 改 
正 后 ，GRACE-C 和 GRACE-D 卫星 相位 残 差 最 大 值 分 别 减 小 至 2.9 mm 和 2.8 mm; 在 高 
度 角 309 ~ 909 范围 内 ， 载 波 相 位 残 差 小 于 1.5 mm。 载 波 残 差 体现 了 误差 模型 和 所 使 用 数 
据 的 一 致 性 ， 可 以 推断 天 线 PCV 的 改正 将 有 助 于 提升 定 轨 精度 。 
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YE: a) GRACE-C 卫星 ; b) GRACE-D 卫星 。 
2 Ke PCV 改正 后 GRACE-FO 卫星 载波 相位 残 差 
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4.2 ”模糊 度 国定 效果 

利用 取 整 法 对 浮 点 模糊 度 进行 固定 ， 宽 起 、 罕 埠 模 糊 度 固 定 阔 值 分 别 为 0.25 和 0.15 
周 。 对 利用 CODE 产品 解 算 的 宽 巷 和 窗 埠 模糊 度 小 数 部 分 进行 统计 ， 如 图 3 所 示 ， 宽 埠 模 
糊 度 小 数 部 分 均值 为 0.01 周 ， 标 准 差 为 0.11 Al; 窗 埠 模糊 度 均 值 为 0， 标准 差 为 0.06 周 。 
较 小 的 标准 差 有 利于 模糊 度 固 定 ， 也 反映 了 相位 偏差 产品 质量 较 高 以 及 所 使 用 模型 的 正确 
性 。 基 于 CNES fü WHU 产品 的 宽 埠 、 罕 起 模 糊 度 小 数 部 分 结果 类 似 。 


8.0 = | 8.0r | = " 
有 均值 一 0.01 标准 差 =0.11 | 有 四 均值 二 0.00 fiif 2:—0.06 | 
6.0; 
o BS 
s : 
S 40 ZR 
xs um 
2.0 
0.0 : 
-0.4 -0.2 0 0.2 0.4 -0.4 -0.2 0 0.2 0.4 
宽 埠 模糊 度 小 数 部 分 / 周 宕 巷 模 糊 度 小 数 部 分 / 周 
a) b) 


注 : a) 宽 埠 模糊 度 小 数 部 分 ，b) 窄 埠 模 糊 度 小 数 部 分 。 
3 ”利用 CODE 产品 解 算 的 GRACE-FO 卫星 模糊 度 小 数 部 分 分 布 


GRACE-FO 卫星 宽 埠 和 和 罕 埠 模糊 度 的 固定 率 如 图 4 Pitas, a) fil c) 为 GRACE-C TE 
宽 巷 、 罕 埠 模 糊 度 固定 率 ，b) Ald) 为 GRACE-D 卫星 宽 埠 、 窄 埠 模 糊 度 固定 率 。 红 色 、 绿 
色 和 蓝 色 点 分 别 表示 使 用 CODE, CNES 和 WHU 产品 的 结果 ， 横 坐标 为 年 积 日 。 年 积 
33 至 51 日 ， 由 于 数据 缺失 未 给 出 结果 。 可 以 看 出 ， 基 于 不 同 产品 的 宽 玲 模糊 度 固定 率 均 优 
于 9996, ^UE REEF 9596, EH CODE 产品 的 模糊 度 固 定 率 最 高 ， 宽 埠 、 罕 埠 
模糊 度 固定 率 分 别 优 于 99.496 和 96.896. 
4.3 ”绝对 定 轨 结 果 分 析 

GRACE-FO 数据 包 中 同时 提供 了 JPL 解 算 的 精密 轨道 产品 ， 该 轨道 是 基于 单 星 模 糊 
度 固 定 的 简化 动力 学 轨道 。 将 文中 生成 的 定 轨 结 果 与 IPL 精密 轨道 进行 比较 ， 以 评估 
定 轨 精 度 。 首 先 ， 为 了 验证 天 线 PCV 的 改正 效果 ， 分 别 将 采用 和 不 采用 PCV 改正 的 定 
轨 结 果 与 JPL 精密 轨道 作 差 。 图 5 为 天 线 PCV 改正 前 后 利用 CODE 产品 的 定 轨 结 果 与 
JPL 轨道 切 向 差异 的 均 方 根 误差 (root mean square, RMS) 值 。a) 和 b) 分 别 为 GRACE-C 
和 GRACE-D 卫星 的 结果 ， 蓝 色 点 为 无 PCV 改正 的 差 值 ， 红 色 点 表示 采用 了 PCV KE. 
经 天 线 PCV 改正 后 ， 两 颗 卫 星 的 切 向 差 值 分 别 由 7.2 mm 和 6.4 mm 减 小 至 6.0 mm 和 
5.9 mm， 径 向 和 法 向 结果 见 表 2， 天 线 PCV 改正 对 轨道 切 向 精度 提升 最 为 显著 。 
图 6 分 别 为 基于 不 同 产品 的 模糊 度 固定 解 轨道 与 JPL 精密 轨道 的 切 向 差 值 ，a) 和 b) 
分 别 对 应 GRACE-C 和 GRACE-D 了 卫星， 红色、 绿色 和 蓝 色 点 分 别 表示 利用 CODE, CNES 
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x xX 10 
a a 
13 R 9 
| m | 
4a 80 “CODE=99.5% -CNES=99.4% “WHU=99.4% | z 80 *CODE=99.4% “CNES=99.4% “WHU=99.3% 
$ 80 100 200 300 400 0 100 200 300 400 
年 积 日 /d 年 积 日 /d 
a) b) 
sc OD dept tia f et ur CURES tap ae ee üt 
B Be: SRN ALES m ES a ER ET E EE T ERE 
H a oa tay E pel ci , R 90 ES "ai CES Nu 2G EE E LM 
iil "ues we | BH | NI 
i 80 -CODE= 96.8% »CNES=95.1% -WHU=96.5% n 80 *CODE- 97.5% : CNES-— 95.8% "WHU- 97.2% 
0 100 200 300 400 0 100 200 300 400 
年 积 日 /d 年 积 日 /d 
c) d) 
ik: a) GRACE-C 卫星 宽 埠 固定 率 ; b) GRACE-D BASES EX; c) GRACE-C 卫星 窗 埠 固定 率 ; d) 
GRACE-D IDEs. 
图 4 GRACE-FO 卫星 模糊 度 固 定 率 
A 30 | ! — A 30 
— “PCV 改 正 =6.00 .无 PCV 改 正 =7.18 S | “PCV 改 正 =5.91 .无 PCV 改 正 =6.39 
WN o 
= 20r 3 . - E: 20 . 
Rd | : ed f 
Es E "v. eet ose PG NECS bz 10 x Tm "tls H £r Y aa S pem TE 
= ry. weigh OA ARE a e Ear NNI IUE YY 
E 0 100 200 300 400 = 0 100 200 300 400 
= 年 积 日 /d * 年 积 日 /d 
al) b) 
ik: a) GRACE-C FÆ; b) GRACE-D 卫星 。 
图 5 PCV 改正 前 后 固定 解 轨道 与 JPL 精密 轨道 的 切 向 差 值 
#2 GRACE-FO 天 线 PCV 改正 对 定 轨 结 果 的 影响 
卫星 轨道 分 量 ”无 PCV 固定 解 /mm PCV 改正 固定 解 /mm PCV 改正 提升 效果 /% 
转向 7.2 6.0 16.7 
GRACE-C 法 向 6.2 6.2 0 
径 向 5.2 5.1 1.9 
转向 6.4 5.9 7.8 
GRACE-D 法 向 6.3 6.1 3.2 
径 向 5.1 5.1 0 


和 WHU 产品 的 结果 。 定 轨 过 程 中 均 利用 4.1 节 结 果 对 天 线 PCV 进行 了 改正 。 法 向 和 径 向 
统计 结果 见 表 3。 可 以 看 出 ， 基 于 不 同 产 品 的 GRACE-FO 模糊 度 固 定 解 轨 道 与 JPL 精密 
轨道 的 差 值 均 小 于 7.0 mm， 不 同 机 构 产 品 定 轨 结果 差异 小 于 1.0 mm。 图 中 部 分 时 间 如 年 积 
H 103, 134 和 270 d 切 向 定 轨 误差 较 大 ， 主 要 是 由 于 该 时 段 内 可 视 GPS 卫星 数 较 少 ， 导 致 
定 轨 精 度 变 差 。 

GRACE-FO 卫星 搭载 了 SLR 激光 角 反 射 器 ，SLR 测 距 值 可 作为 独立 手段 对 卫星 定 


ive 
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30 r r | 30, r = r 
g -CODE-6.0 -CNES-6.7 : WHU=6.7 g - CODE-5.9 -CNES—6.5 - WHU—6.6 
A 20| = | 20} | 
= 10 S rol 
Ey sti ARS AER ON X RV salon RAAN ES 
0 100 200 300 400 Qj 100 200 300 400 
年 积 日 /d 年 PA /d 
a) b) 
iE: a) GRACE-C 卫星 ， b) GRACE-D 卫星 。 
6 ”基于 不 同 机 构 产 品 的 固定 解 轨 道 与 JPL 精密 轨道 的 切 向 差 值 
表 3 基于 不 同 机 构 产 品 的 固定 解 轨道 与 JPL 精密 轨道 的 三 轴 差 值 
卫星 产品 源 ” 切 向 /mm 法 向 /mm 径 向 /mm _ 三 轴 平 均 /mm 
CODE 6.0 6.2 5.1 5.8 
GRACE-C CNES 6.7 6.4 5.0 6.1 
WHU 6.7 6.5 5.3 6.2 
CODE 5.9 6.1 5.1 5.7 
GRACE-D CNES 6.5 6.3 4.9 5.9 
WHU 6.6 6.5 5.3 6.2 
轨 结 果 进 行 评 佑 。 选 取 数 据 质量 较 好 的 9 SLR 测 站 (Yarragadee, Greenbelt, Haleakala, 


Zimmerwald, Mount Stromlo, Wettzell, Graz, Herstmonceux, Potsdam) 进行 分 析 。 利 用 解 
算 的 GRACE-FO 卫星 轨道 和 SLR 测 站 坐标 计算 测 站 与 卫星 间 的 距离 ， 并 与 SLR 观测 值 作 


差 ， 其 残 差 可 反映 定 轨 精 度 。 图 7 为 SLR 测 距 残 差 值 ， 大 于 20 cm RAE AA al BR 
图 7a), b) 分 别 表 示 GRACE-C 和 GRACE-D 卫星 结果 ， 


红色 、 绿 色 和 蓝 色 点 分 别 表示 利 
CODE、CNES 和 WHU 产品 的 结果 。 利 用 CODE 和 WHU 产品 的 定 轨 结 果 与 SLR FF 
合 较 好 ， 检 核 精 度 分 别 优 于 9.6 mm 和 9.1 mm， 基 于 CNES 产品 的 检 核 精度 优 于 10.7 mm, 
SLR 检验 精度 整体 差异 小 于 2 mm. 


引 200 = 200, 
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a) b) 
YE: a) GRACE-C 卫星 ; b) GRACE-D 卫星 。 
图 7 基于 不 同 机 构 产 品 的 固定 解 轨 道 SLR 检 核 残 差 
4.4 星 间 基线 结果 分 析 
GRACE-FO 卫星 同时 搭载 了 KK 波段 激光 测 距 仪 ， 其 测 距 精度 可 达 微 米 级 号。 利用 单 
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精密 定 轨 结 果 作 差 计 算出 两 颗 卫 星 间 的 距离 ， 与 KBR 观测 值 进 和 
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道 的 切 向 精度 。 图 8 为 2019 年 180 d 利用 精密 轨道 计算 的 星 间 距离 与 KBR 测 距 值 差 值 的 
标准 差 (standard deviation, STD)， 其 中 红色 点 为 采用 PCV 改正 的 浮 点 解 结果 ， 绿 色 为 无 
PCV 改正 的 固定 解 结果 ， 蓝 色 为 采用 PCV 改正 的 固定 解 结果 。 可 以 看 出 ， 固 定 模糊 度 可 
以 将 测 距 残 差 由 5.4 mm 减 小 至 1.7 mm， 对 天 线 PCV 进行 改正 可 以 将 测 距 残 差 由 3.0 mm 
减 小 至 1.7 mm。 图 中 浮 点 解 结果 在 12 h 左右 误差 较 大 ， 主 要 是 由 于 该 时 段 可 视 GPS 卫星 
数 较 少 ， 正 确 固 定 模糊 度 能 够 有 效 提升 该 时 段 的 定 轨 精度 。 图 9 为 2019 年 全 年 KBR 测 距 
残 差 的 STD 值 。 利 用 CODE, CNES 和 WHU 产品 解 算 的 GRACE-FO 卫星 精密 轨道 星 间 
KBR 残 差 STD 值 分 别 优 于 1.8, 2.3 和 2.1 mm， 与 文献 [26] 中 评估 结果 类 似 。 
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8 ”利用 2019 年 180 d 轨道 计算 的 KRB 测 距 9 ”基于 不 同 机 构 产 品 的 固定 解 轨道 KBR RE 
RE 标准 差 


5 结 论 


本 文 利 用 宽 埠 和 窄 埠 相位 偏差 产品 对 载波 相位 观测 值 进行 改正 ， 实 现 了 单 星 精密 定 轨 
中 模糊 度 参数 的 固定 ， 同 时 基于 无 电离 层 组 合 载波 相位 残 差 估 计 天 线 PCV 改正 数 ， 实 现 了 
星 载 接收 天 线 的 在 轨 标 定 。 利 用 GRACE-FO 卫星 2019 年 在 轨 实 测 数据 ， 分 别 采 用 CODE, 
CNES 和 WHU 发 布 的 GPS 精密 轨道 、 钟 差 和 相位 偏差 产品 进行 GRACE-FO 卫星 模糊 
度 固定 和 精密 定 轨 。 结 果 表 明 ， 天 线 PCV 在 轨 标 定 能 够 有 效 减 小 载波 相位 残 差 ， 可 进 
步 提升 定 轨 精度 ， 其 中 对 轨道 切 癌 提升 最 为 明显 。 采 用 不 同 机 构 产 品 的 GRACE-FO 卫星 
宽 埠 模糊 度 固定 率 优 于 99 允 ， 鹤 巷 模 糊 度 固定 率 优 于 95%， 其 中 基于 CODE 产品 的 模糊 度 
固定 率 最 高 ， 宽 埠 固 定 率 优 于 99.4%, FART 96.8% DÀ JPL 精密 轨道 作为 参考 ， 利 用 
不 同 机 构 产 品 解 算 的 GRACE-FO 精密 轨道 三 轴 平 均 精 度 优 于 7.0 mm。 利 用 SLR 测 距 值 
对 GRACE-FO 定 轨 结 果 进 行 检验 ， 其 中 基于 CODE 和 WHU 产品 的 定 轨 结果 SLR RH 
RMS 值 小 于 10.0 mm, CNES 产品 定 轨 结 果 SLR 残 差 优 于 11.0 mm. S HER] KBR W 
距 值 对 GRACE-FO 卫星 相对 轨道 精度 进行 评估 ， 基 于 CODE, CNES 和 WHU 产品 的 相对 
轨道 精度 分 别 优 于 1.8, 2.3 812.1 mm。 不 同 机 构 产 品 GRACE-FO 卫星 精密 定 轨 结果 整体 
差异 小 于 2 mm。 另外 ， 在 本 文 分 析 中 ， 部 分 观测 时 段 内 由 于 可 视 卫 星 数 较 少 导致 基于 GPS 
的 简化 动力 学 定 轨 精 度 变 差 ， 后 续 随 着 多 模 GNSS 星 载 接收 机 陆续 在 轨 应 用 ， 结 合 多 系统 
相位 偏差 产品 ， 基 于 模糊 度 固 定 的 单 星 定 轨 精 度 将 进一步 提升 。 
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GRACE-FO Precise Orbit Determination with Single 


Receiver Ambiguity Resolution 


JIN Biao!?2?, LI Yu-qiang!, ZHOU Weit, LI Zhu-lian', CHEN Shan-shan? 


(1. Yunnan Observatories, Chinese Academy of Sciences, Kunming 650216, China; 2. University 
of Chinese Academy of Sciences, Beijing 100049, China; 3. Space Star Technology Co., Ltd., Beijing 
100094, China; 4. Beijing Institute of Tracking and Telecommunications Technology, Beijing 100094, 
China; 5. Beijing Sixents Technology Co., Ltd., Beijing 100094, China) 


Abstract: Single receiver ambiguity resolution can significantly improve the accuracy of 
orbit determination for low earth orbit (LEO) satellites. Carrier phase bias product is em- 
ployed to calibrate the satellite phase delay, and single difference (SD) observation between 
satellites is applied to remove the phase delay of receiver. Combined with related GPS pre- 
cision orbit and clock products, wide lane and narrow lane ambiguities are fixed to integer. 
Ionosphere free (IF) ambiguity is recovered with the fixed wide lane and narrow lane am- 
biguities. The resolved SD ionosphere free ambiguities are taken as pseudo observations to 
constrain the undifferenced IF ambiguity parameters and derive the orbit of LEO satellite. 
Antenna phase center variation (PCV) is estimated with carrier phase residuals of IF com- 
bination. GPS precise orbits, clock and phase bias products provided by Center for Orbit 
Determination in Europe (CODE), Centre National d'Etudes Spatiales (CNES) and Wuhan 
University (WHU) are used to perform the single receiver ambiguity resolution and precise 
orbit determination. The performance of ambiguity resolution and orbit accuracy related to 
the different products is analyzed. In-flight antenna PCV correction can effectively reduce 
the carrier phase residuals and improve the orbit accuracy, especially the along-track orbit 
accuracy. With products from different institutions, the fixing rates of wide-lane and narrow- 
lane ambiguities are better than 9996 and 9596. Compared with the precise orbit generated 
by Jet Propulsion Laboratory (JPL), the one dimensional accuracy of GRACE-FO's orbit 
determined with single receiver ambiguity resolution is better than 7.0 mm. RMSs of SLR 
residuals for GRACE-FO orbits generated with CODE, CNES and WHU's products are 9.6, 
10.7 and 9.1 mm. Standard deviations (STDs) of the independent K-band range (KBR) 
residuals for the three kinds of orbits are 1.8, 2.3 and 2.1 mm respectively. The overall 
difference of GRACE-FO orbits calculated with different products is less than 2 mm. 


Key words: GRACE-FO satellites; single receiver ambiguity resolution; antenna phase 


center variation (PCV); precise orbit determination 


